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Abstract 

In a balloon drawing of a tree, all the children under the same parent are 
placed on the circumference of the circle centered at their parent, and the 
radius of the circle centered at each node along any path from the root 
reflects the number of descendants associated with the node. Among various 
styles of tree drawings reported in the literature, the balloon drawing enjoys 
a desirable feature of displaying tree structures in a rather balanced fashion. 
For each internal node in a balloon drawing, the ray from the node to each 
of its children divides the wedge accommodating the subtree rooted at the 
child into two sub-wedges. Depending on whether the two sub-wedge angles 
are required to be identical or not, a balloon drawing can further be divided 
into two types: even sub-wedge and uneven sub-wedge types. In the most 
general case, for any internal node in the tree there are two dimensions of 
freedom that affect the quality of a balloon drawing: (1) altering the order in 
which the children of the node appear in the drawing, and (2) for the subtree 
rooted at each child of the node, flipping the two sub-wedges of the subtree. 
In this paper, we give a comprehensive complexity analysis for optimizing 
balloon drawings of rooted trees with respect to angular resolution, aspect 
ratio and standard deviation of angles under various drawing cases depending 
on whether the tree is of even or uneven sub- wedge type and whether (1) and 
(2) above are allowed. It turns out that some are NP-complete while others 
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can be solved in polynomial time. We also derive approximation algorithms 
for those that are intractable in general. 

Key words: tree drawing, graph drawing, graph algorithms 



1. Introduction 

Graph drawing addresses the issue of constructing geometric represen- 
tations of graphs in a way to gain better understanding and insights into 
the graph structures. Surveys on graph drawing can be found in [il, 3. If 
the given data is hierarchical (such as a file system), then it can often be 
expressed as a rooted tree. Among existing algorithms in the literature for 
drawing rooted trees, the work of 
ing binary trees. The idea behind 



11] developed a popular method for draw- 
ll| is to recursively draw the left and right 
subtrees independently in a bottom-up manner, then shift the two drawings 
along the x-direction as close to each other as possible while centering the 
parent of the two subtrees one level up between their roots. Different from 



tne conventional triangular tree arawmg 01 ftu-arawmgs [iz|j, raaiai 
drawings jsl and balloon drawings j2l, 0,13, lol. Iiof are also popular for vi- 
sualizing hierarchical graphs. Since the majority of algorithms for drawing 
rooted trees take hnear time, rooted tree structures are suited to be used in 
an environment in which real-time interactions with users are frequent. 

Consider Figured] for an example. A balloon drawing ji], 0, [oj of a rooted 
tree is a drawing having the following properties: 

• all the children under the same parent are placed on the circumference 
of the circle centered at their parent; 

• there exist no edge crossings in the drawing; 

• the radius of the circle centered at each node along any path from the 
root node reflects the number of descendants associated with the node 
(i.e., for any two edges on a path from the root node, the farther from 
the root an edge is, the shorter its drawing length becomes). 

In the balloon drawing of a tree, each subtree resides in a wedge whose 
end-point is the parent node of the root of the subtree. The ray from the 
parent node to the root of the subtree divides the wedge into two sub-wedges. 
Depending on whether the two sub-wedge angles are required to be identical 
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Figure 1: Illustration of balloon drawings with (a) even sub- wedges and (b) uneven sub- 
wedges, where each node is drawn by a point; each edge is drawn by a solid straight line 
segment; the center of the largest circle in a balloon drawing is the root node. 

or not, a balloon drawing can further be divided into two types: drawings 
with even sub-wedges (see Figure [U^a)) and drawings with uneven sub-wedges 
(see Figure dl^b)). One can see from the transformation from Figure [T]^a) to 
Figure [U^b) that a balloon drawing with uneven sub- wedges is derived from 
that with even sub-wedges by shrinking the drawing circles in a bottom-up 
fashion so that the drawing area is as small as possible joj. Another way to 
differentiate the two is that for the even sub-wedge case, it is required that 
the position of the root of a subtree coincides with the center of the enclosing 
circle of the subtree. 

Aesthetic criteria specify graphic structures and properties of drawing, 
such as minimizing number of edge crossings or bends, minimizing area, 
and so on, but the problem of simultaneously optimizing those criteria is, 
in many cases, NP-hard. The main aesthetic criteria on the angle sizes in 
balloon drawings are angular resolution, aspect ratio, and standard deviation 
of angles. Note that this paper mainly concerns the angle sizes, while it is 
interesting to investigate other aesthetic criteria, such as the drawing area, 
total edge length, etc. Given a drawing of tree T, an angle formed by the two 
adjacent edges incident to a common node v is called an angle incident to 
node V. Note that an angle in a balloon drawing consists of two sub- wedges 
which belong to two different subtrees, respectively (see Figure [T]). With 
respect to a node v, the angular resolution is the smallest angle incident to 
node V, the aspect ratio is the ratio of the largest angle to the smallest angle 
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incident to node v, and the standard deviation of angles is a statistic used as 
a measure of the dispersion or variation in the distribution of angles, equal 
to the square root of the arithmetic mean of the squares of the deviations 
from the arithmetic mean. 

The angular resolution (resp., aspect ratio; standard deviation of angles) 
of a drawing of T is defined as the minimum angular resolution (resp., the 
maximum aspect ratio; the maximum standard deviation of angles) among 
all nodes in T. The angular resolution (resp., aspect ratio; standard deviation 
of angles) of a tree drawing is in the range of (0°,360°) (resp., [l,C)o) and 
[0, oo)). A tree layout with a large angular resolution can easily be identified 
by eyes, while a tree layout with a small aspect ratio or standard deviation 
of angles often enjoys a very balanced view of tree drawing. It is worthy of 
pointing out the fundamental difference between aspect ratio and standard 
deviation. The aspect ratio only concerns the deviation between the largest 
and the smallest angles in the drawing, while the standard deviation deals 
with the deviation of all the angles. 

With respect to a balloon drawing of a rooted tree, changing the order 
in which the children of a node are listed or flipping the two sub-wedges of 
a subtree affects the quality of the drawing. For example, in comparison 
between the two balloon drawings of a tree under different tree orderings 
respectively shown in Figures [2]^a) and[2]^b), we observe that the drawing in 
Figure MJo) displays little variations of angles, which give a very balanced 
drawing. Hence some interesting questions arise: How to change the tree 
ordering or flip the two sub-wedge angles of each subtree such that the balloon 
drawing of the tree has the maximum angular resolution, the minimum aspect 
ratio, and the minimum standard deviation of angles ? 

Throughout the rest of this paper, we let RE, RA, and DE denote the 
problems of optimizing angular resolution, aspect ratio, and standard devi- 
ation of angles, respectively. In this paper, we investigate the tractability of 
the RE, RA, and DE problems in a variety of cases, and our main results 
are listed in Table [H in which trees with 'flexible' (resp., 'fixed') uneven 
sub- wedges refer to the case when sub- wedges of subtrees are (resp., are not) 
allowed to flip; a 'semi-ordered' tree is an unordered tree where only the 
circular ordering of the children of each node is fixed, without specifying if 
this ordering is clockwise or counterclockwise in the drawing. Note that a 
semi-ordered tree allows to flip uneven sub-wedges in the drawing, because 
flipping sub-wedges of a node in the bottom-up fashion of the tree does not 
modify the circular ordering of its children. See Figure [2] for an experimental 
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(c) Initail drawing with uneven sub-wedges (d) RA2 




(e) RA3 (f) RA4 



Figure 2: An experimental example, where (a) and (c) are initial balloon drawings with 
even and uneven sub- wedges, respectively; (b), (d), (e) and (f) achieve the optimality of 
RAl, RA2, RA3 and RA4, respectively. Note that the differences of (d) from (c) are 
encompassed by shaded regions. 
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Table 1: The time complexity for optimizing main aesthetic criteria of balloon drawing. 



aesthetic criterion denotation complexity reference 



CI: 



unordered trees with 

even sub-wedges 



angular resolution REl 0(n log n) 

aspect ratio RAl 0(n log n) 

standard deviation DEI 0{n log n)* 



m 

[Thmg 
[Thm|2] 
[Thm|3] 
[Thmg 
[Thm[5] 
[Thm[6l[8] 
[ThmEigol 

[Thm[5] 
[Thm[6l[8] 
[ThmElllOl 



C2: 



semi-ordered trees with 
flexible uneven sub-wedges 



angular resolution RE2 
aspect ratio RA2 
standard deviation DE2 



Oin)* 
0(n2)* 



C3: 



unordered trees with 
fixed uneven sub-wedges 



angular resolution RES O(nlogn)* 

aspect ratio RA3 NPC* 

standard deviation DE3 NPC* 



C4: 



unordered trees with 
flexible uneven sub-wedges 



angular resolution RE4 O(rilogn)* 

aspect ratio RA4 NPC* 

standard deviation DE4 NPC* 



*Thc marked entries are the contributions of this paper. 
Q for RE2 and RA2 require 0(n^-^) time. 



Note that earlier results reported in 



example with the drawings which achieve the optimahty of RA1-RA4. In 
Table [Tj with the exception of REl and RAl (which were previously ob- 
tained by Lin and Yen in j^), all the remaining results are new. We also give 
2-approximation algorithms for RA3 and RA4, and 0(\/n)-approximation 
algorithms for DE3 and DE4. Finding improved approximation bounds for 
those intractable problems remains an interesting open question. 

The rest of the paper is organized as follows. Some preliminaries are 
given in Section [2l The problems for cases CI and C2 are investigated in 
Section [31 The problems for cases C3 and C4 are investigated in Section HJ 
The approximation algorithms for those intractable problems are given in 
Section [5l Finally, a conclusion is given in Section [6l 

2. Preliminaries 

In this section, we first introduce two conventional models of balloon 
drawing, then define our concerned problems, and finally introduce some 
related problems. 

2.1. Two Models of Balloon Drawing 

There exist two models in the literature for generating balloon drawings 
of trees. Given a node v, let r{v) be the radius of the drawing circle centered 
at V. If we require that r{v) = r{w) for arbitrary two nodes v and w that are 
of the same depth from the root of the tree, then such a drawing is called 
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a balloon drawing under the fractal model [7|]. The fractal drawing of a tree 
structure means that if and r^-i are the lengths of edges at depths m 
and m — 1, respectively, then = 7 x Vm-i where 7 is the predefined ratio 
(0 < 7 < 1) associated with the drawing under the fractal model. Clearly, 
edges at the same depth have the same length in a fractal drawing. 

Unlike the fractal model, the subtrees with nonuniform sizes (abbreviated 
as SNS) model 0,i allows subtrees associated with the same parent to reside 
in circles of different sizes (see also Figure [U^a)), and hence the drawing based 
on this model often results in a clearer display on large subtrees than that un- 
der the fractal model. Given a rooted ordered tree T with n nodes, a balloon 
drawing under the SNS model can be obtained in 0{n) time (see [2I, 0]) 
bottom-up fashion by computing the edge length r and the angle 6*, between 
two adjacent edges respectively according to r = C/(27r) = (2 -Ri)/(27r) 
and 9i = {Ri + freeware + Ri^i)/r (see Figure [^a)) where r is the radius 
of the inner circle centered at node cq; C is the circumference of the inner 
circle; Ri is the radius of the outer circle enclosing all subtrees of the i-th 
child of Co, and Rq is the radius of the outer circle enclosing all subtrees of 
Co; since there exists a gap between C and the sum of all diameters, we can 
distribute to every 6i the gap between them evenly, which is called a free arc, 
denoted by freeware. 

Note that the balloon drawing under the SNS model is our so-called 
balloon drawing with even sub-wedges. A careful examination reveals that 
the area of a balloon drawing with even sub-wedges (generated by the SNS 
model) may be reduced by shrinking the free arc between each pair of subtrees 
and shortening the radius of each inner circle in a bottom- up fashion by 
which we can obtain a smaller-area balloon drawing with uneven sub-wedges 
(e.g., see the transformation from Figure [D^a) to Figure [T|^c)). 

2.2. Notation and Problem Definition 

In what follows, we introduce some notation, used in the rest of this 
paper. A circular permutation vr is expressed as: vr = (tti, 7r2, 7r„) where 
for i = 1,2, ...,n, iTi is placed along a circle in a counterclockwise direction. 
Note that 7r„ is adjacent to tti; i © 1 denotes i + 1 {mod n); i Q 1 denotes 
i — 1 {mod n). Due to the hierarchical nature of trees and the ways the 
aesthetic criteria (measures) for balloon drawings are defined, an algorithm 
optimizing a star graph can be applied repeatedly to a general tree in a 
bottom-up fashion , yielding an optimum solution with respect to a given 
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Figure 3: Notations used in a balloon drawing of a star graph with uneven sub-wedges. 

aesthetic criterion. Thus, it suffices to consider the balloon drawing of a star 
graph when we discuss these problems. 

A star graph is characterized by a root node Cq together with its n children 
Ci, ...,Cn, each of which is the root of a subtree located entirely in a wedge, 
as shown in Figure [It^a) (for the even sub-wedge type) and Figure [3] (for the 
uneven sub- wedge type). In what follows, we can only see Figure [3] because 
the even sub-wedge type can be viewed as a special case of the uneven sub- 
wedge type. The ray from cq to Cj further divides the associated wedge 
into two sub- wedges SWi^ and SWi^i with sizes of angles wo(0 wi{i), 
respectively. Note that wq{i) and wi{i) need not be equal in general. An 
ordering of Cq's children is simply a circular permutation a = {ai, (J2, (Jn), 
in which o"j G {1, 2, n} for each i. 

There are two dimensions of freedom affecting the quality of a balloon 
drawing for a star graph. The first is concerned with the ordering in which 
the children of the root node cq are drawn. With a given ordering, it is also 
possible to alter the order of occurrences of the two sub-wedges associated 
with each child of the root. With respect to child q and its two sub- wedges 
SWifl and SWi^i, we use tj G {0, 1} to denote the index of the first sub-wedge 
encountered in a counterclockwise traversal of the drawing. For convenience, 
we let = 1 — tj. We also write t = (ti,...,tn) (ti G {0, 1}, 1 < i < n), 
which is called the sub-wedge assignment (or simply assignment) . As shown 
in Figure |3l the sequence of sub-wedges encountered along the cycle centered 
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at Co in a counterclockwise direction can be expressed as: 



(1) 



If Wo{i) = Wi{i) for each i G {1, then the drawing is said to be of 

even sub-wedge type; otherwise, it is of uneven sub-wedge type. As mentioned 
earher, the order of the two sub-wedges associated with a child (along the 
counterclockwise direction) affects the quality of a drawing in the uneven sub- 
wedge case. For the case of uneven sub- wedge type, if the assignment t is 
given a priori, then the drawing is said to be of fixed uneven sub-wedge type; 
otherwise, of flexible uneven sub- wedge type (i.e., t is a design parameter). 

As shown in Figure [3], with respect to an ordering a and an assignment 
t in circular permutation ([1]), c^. and Co-^^i, 1 < ^ < "n., are neighboring 
nodes, and the size of the angle formed by the two adjacent edges CqcJ^ and 
coCaiJ^ is 9i = Wt'.{cri) + i^j.^-^ ((jj©i). Hence, the angular resolution (denoted 
by AngResla,t), the aspect ratio (denoted by AspRatio^^t), and the standard 
deviation of angles (denoted by StdDev^^t) can be formulated as 

AngResla,t= min 9i= min {wf ((7^) + Wt^f^i (o-^ei)}; 

l<i<n l<i<n ^ 

. „ ^. maxi<,<„6ii maxi<,<„{u;t;(CTi) +u;t,3i(criei)} 
AspKanOc^t ^ ^ 



mini<i<„6'i mini<i<„{u;4; (ai) + Wt,ei (o'^ei)} 

StdDeVa.t — \l 



'Er=i(«^t^('^0^ + wt.g,i(criei)^) , '^J27=lWt[icr^)wt,^Xa,(sl) f^jA^ 



n n \ n 



We observe that the first and third terms inside the square root of the above 
equation are constants for any circular permutation a and assignment t, and 
hence, the second term inside the square root is the dominant factor as far 
as StdDev^^t is concerned. We denote by SOP^ t the sum of products of 
sub-wedges, which can be expressed as: 



We are now in a position to define the RE, RA and DE problems in Ta- 
ble [T] for four cases (CI, C2, C3, and C4) in a precise manner. The four 
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cases depend on whether the circular permutation a and the assignment t in 
a balloon drawing arc fixed (i.e., given a priori) or ficxiblc (i.e., design pa- 
rameters). For example, case C3 allows an arbitrary ordering of the children 
(i.e., the tree is unordered), but the relative positions of the two sub- wedges 
associated with a child node are fixed (i.e., fiipping is not allowed). The 
remaining three cases are easy to understand. 

We consider the most fiexible case, namely, C4, for which both a and t 
are design parameters, which can be chosen from the set S of all circular per- 
mutations of {1, ...,n} and the set T of all n-bit binary strings, respectively. 
The RE and RA problems, respectively, are concerned with finding a and t 
to achieve the following: 

optAngResl = umx^{AngResl„^t};optAspRatio = imn ^{AspRatiOfj^t} ■ 

The DE problem is concerned with finding a and t to achieve the following: 

optStdDev = min {StdDeVat}- 
o-eS;teT 

As stated earlier, optStdDev is closely related to the SOP problem, which is 
concerned with finding a and t to achieve the following: 

optSOP = min {SOP^^t}- 

2.3. Related Problems 

Before deriving our main results, we first recall two problems, namely, 
the two-station assembly line problem (2SAL) and the cyclic two-station 
workforce leveling problem (2SLW) that are closely related to our problems 
of optimizing balloon drawing under a variety of aesthetic criteria. Con- 
sider a serial assembly line with two stations, say STi and ST21 and a set 
J = { Ji, J2, J„} of n jobs. Each job Jj = {Wn, Wi2) consists of two tasks 
processed by the two stations, respectively, where Wn (resp., Wi2) is the 
workforce requirement at STi (resp., ST2). Assume the processing time of 
each job at each station is the same, say r/n. Consider a circular permuta- 
tion (J^j, J52, J5„) of J where 5 = {61,62, ...,Sn) is a circular permutation 
of {1, 2, n}. At any time point, a single station can only process one job. 
We also assume that the two stations are always busy. During the first time 
range [0,r/n], Jg^ and are processed by ST2 and STi, respectively, and 
the workforce requirement is Ws^2 + Ws^i- Similarly, for each i, during the 
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^ jfr] [sT^ 

— ^ — — 

S J4 rJSSl 



time range 




ST2 


workforce requirement 


[ 0, r/4] 


J2 


J3 


2 + 1=3 


[T/4,2r/4] 


J4 


J2 


7 + 4 = 11 


[2r/4,3T/4] 


Ji 


J4 


2 + 2= 4 


[3r/4,4T/4] 


J3 


Ji 


3 + 6= 9 



Figure 4: An example for 2SAL and 2SLW. 



time range [{i — l)r/?2, ir/ri], J^. and J^.^^ are processed at ST2 and 5'Ti 
stations respectively, and the workforce requirement is Ws^2 + W^^^eii- 

For example, consider J = { Ji, J2, J3, J4} where Ji = (2, 3), J2 = (1, 7), 
J3 = (6, 2), and J4 = (4, 2). For a certain circular permutation (J3, J2, J4, Ji) 
of J, the workforce requirements for each period of time as well as the jobs 
served at the two stations are given in Figure HI where the largest workforce 
requirement is 11; the range of the workforce requirements among all the 
time periods is [3,11]. 

The 25*^4^ and 2SLW problems are defined as follows: 

• 2SAL: Given a set of n jobs, find a circular permutation of the n jobs 
such that the largest workforce requirement is minimized. 

• 2SLW (decision version): Given a set of n jobs and a range [LB, UB] 
of workforce requirements, decide whether a circular permutation exists 
such that the workforce requirement for each time period is between 
LB and UB. 

It is known that 2SAL is solvable in 0(?T,log?T,) time ^], while 2SLW is NP- 



complete [13 



3. Cases CI (Unordered Trees with Even Sub-Wedges) and C2 
(Semi-Ordered Trees with Flexible Uneven Sub- Wedges) 

First of all, we investigate the DEI problem (SOPl problem), i.e., finding 
a balloon drawing optimizing optSOP for case CI (i.e., unordered trees with 
even sub- wedges). In this case, the two sub- wedges associated with a child 
node in a star graph are of the same size. For notational convenience, we 
order the set of wedge angles {wo{i) + wi{i) : i = 1, ■ ■ ■ ,n} (note that in this 
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case WQ{i) = Wi{i) for each i) in ascending order as either 

mi,m2,---,mfc_i,mfc,Mfc,Mfc_i,---,M2,Mi if n = 2A;, or (3) 
mi,m2, ■ ■ ■ ,mA;-i,mfc,mz(i, Mfc, Mfc_i, ■ ■ ■ , M2, Ml if n = 2A; + 1, (4) 

for some k, where rrii (resp., Mj) is the i-th minimum (resp., maximum) 
among all, and mid is the median if the number of elements is odd. Note 
that the size of each angle between two edges in the drawing may be one 
of the forms (m« + mb)/2, (m^ + Mb)/2, {Ma + mf,)/2, or (M„ + Mb)/2 for 
some a, 6 G {1, ■ ■ ■ , n}, and hence, there may exist more than one angle with 
the same value. In what follows, we are able to solve the DEI problem by 
applying Procedure [H 

Procedure 1 OptBalloonDrawing-DEI 

Input: a star graph S with n child nodes of nonuniform sizes 

Output: a balloon drawing of S optimizing standard deviation of angles 

1: sort {wo{i) + wi{i) : i = 1, • • • , n} in ascending order as either Equation ([3]), if 

n = 2k, or Equation (jH), if n = 2/c + 1 
2: for convenience, let the child node with wedge m^, mid or Mj be also denoted 

by rrii, mid or Mj, respectively 
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ii n = 2k then 
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if k is odd then 




















5 
6 


output (Ml, m2 
else 


M3 


m4, • 


• , Mfe_i 


,mk,Mk,mk- 




,M4 




M2, 


mi) 


7 

8 


output (Ml, m2 
end if 


M3 


m4, • 


• ,mk-i 
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1) • • 
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else 
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if k is odd then 
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12 


output (Ml, ?n,2, 
else 
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m-4, • • 


■ , Mk-i 


mk,mid, Mk 


nik- 


-ir ■ ■ 


,M4 




,M2,mi 


13 
14 


output (Ml, m2, 
end if 




7714, • • 


■ ,mk-i, 


Mk,mid, mk 


Mk- 


-I7 • • • 


,M4 




,M2,mi 


15 


end if 





















Theorem 1. The DEI problem is solvable in 0{nlogn) time. 

Proof. In what follows, we show that Procedure [H which clearly runs in 
0{nlogn) time, can be applied to correctly producing the optimum solution. 
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We only consider an output case in Procedure [D 



a = (Ml, 1712, M3, rrii, Mfc_i, m^, mid, M^, rrik-i, M4, m^, M2, mi) 

i.e., n = 2/c + l and k is odd; the remaining cases are similar (in fact, simpler). 
Note that SOP^^t = (^\Zi Mimi+i + mk x mid + mid x Mk + Yli=i ^i^i+i + 
miMi)/4, for this output case. 

We proceed by induction on an integer number i, for i = 1 to k, to prove 
that, with respect to the SOP measure, no circular permutations perform 
better than a certain circular permutation S which contains the sequence 



If the above holds, then no circular permutations perform better than a cer- 
tain circular permutation 6 which contains sequence Sk- That is, no circular 
permutations perform better than circular permutation 6 = {Sk,iTLid) = a, 
as required. 

For i = 1, we show that no circular permutations perform better than a 
certain circular permutation 6 which contains sequence Si = mi Mi. Con- 
trarily suppose that there exists a circular permutation 6' in which mi is not 
adjacent to Mi so that SOP51 < SOP5. We assume that mi (resp.. Mi) is 
adjacent to x = mi + li (resp., y = mi + I2) in S' where mi < x,y < Mi, 
X y, and /i,/2 > 0. W.l.o.g., let 6' be {xmiS'yMiS") where S' U S" = 
{m2, ■ ■ ■ , m„, mid, Mn, ■ ■ ■ , M2} \ {x, y}. Consider circular permutation S = 
{xyS'^miMiS") where 5"^ is the reverse of S' . Then SOPs' — SOPs = 
{xmi + yMi - xy - miMi)/4 = kiMi - mi- /i)/4 = kiMi - x)/4 > 0, 
which is a contradiction. 

Suppose that no circular permutations perform better than a certain cir- 
cular permutation which contains sequence Si-i. We show that no circular 
permutation perform better than a certain circular permutation 6i which 
contains sequence S^. In the following, we only consider the case when i is 
even (i.e.. Si = MiSi^imi); the other case is similar. 

Contrarily suppose that there exists a circular permutation 5- which per- 
form better than 6i, i.e., SOPs'_ < SOPs^. By the inductive hypothesis, 
SOPs'_ > SOPs._^ for some circular permutation 6i-i which contains se- 
quence Si^i. W.l.o.g., suppose that 5j_i = {Si^iXiS'miX2S"x^MiS"'x/i) 
where < xi,---,Xi < Mi and S' U S" U S'" = {mi^i, ■ ■ ■ ,mn,mid. 




miMi, if z = 1; 
MiSi-imi, if i is even; 
miSi-iMi, if i is odd. 
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M„, ■ ■ ■ , Mj+i} \ {a;i, ■ ■ ■ , X4}; the other cases are similar. Assume Xi = 
rrii + /i, ■ ■ X4 = mj + 1^ where /i, • • ■ , /4 > 0. Let Mi = mi + where > Ij 
for each j G {1, ■ ■ ■ ,4}. Consider 6i = {Si-imiS'^xiX2S"xzXiS"'^Mi) . Then 
SOP5^_-^-SOPs^ = {Mi_iXi + miX2 + X3Mi+X4mi_i-Mi_imi-xiX2-X3X4- 
Mimi.i)/A = h{Mi.i - - h)/^ + (m,_i - rm - h){h - k)/4: = h{Mi_i - 
a;2)/4+(mi_i-a;3)(/4-/5)/4 > 0. Hence, SOP^^^ > SOPs^_, > SOPs^, which 
is a contradiction. □ 

Now consider case C2 (semi-ordered trees with flexible uneven angles). In 
this case, the ordering of children of the root, cr = (1, 2, ■ ■ ■ , n), is fixed, and 
only the assignment of t = (ti, ■ ■ ■ , t„) needs to be specified. Our solutions for 
RE2, RA2 and DE2 are based on dynamic programming approaches. Those 
results are given as follows: 

Theorem 2. The RE2 problem can he solved in 0{n) time. 

Proof. W.l.o.g., assume a = (1, 2, n). Recall from Equation ([1]) that ift = 
(ti, ...,tn) is the assignment of sub- wedges, then the sequence of sub- wedges 
encountered in a counterclockwise direction is {wt^ (1) , Wt'^{l) , Wt^i^) , Wt'^{2) , ■ ■ 
Wt^{n),wt>^{n)). We define fi{wt^{l),Wt'^{{)) as follows: 



max {mm{{wt: {l)+wt^{2)), {wt'^{2)+wt^{3)), {wf^ {i-l)+WtXi))}}- 

1},2<J<«-1 ^ ^ ' ^ 



That is, the solution maximizes the minimum sum of adjacent sub- wedge 
pairs for the first i children, given Wii(l) and Wt'^{i) as the outer sub-wedges 
of first child and i-th child, respectively. Notice that Wt'.{i) +wi{l) is not in- 
cluded in calculating fi{wt^{l),Wt'^{i)), meaning that the first child is not con- 
sidered to be adjacent to the i-th child. We can observe that fi{wt^{l), Wt'.{i)) 
can be formulated as the following dynamic programming formula: 

f^{wt,{l),wt'{i)) = max {mm{fi_i{wt,{l),wt'^ (i-l)), Wt'^ {i-l)+WtXi)}}- 
Finally, we have: 

optAngResl= max {min{/„(wij(l),t(7t;Jn)), Wt^{l) + Wt'^^{n)}}. 

It is easy to see that the above algorithm gives the correct answer and runs 
in linear time. □ 
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Theorem 3. The RA2 problem can he solved in 0{rP) time. 

Proof. Since only flipping subwedges is allowed in this case, Wo{i) and Wi{i) 
can be the neighbors of Wo{iQ)l) and Wi{i(Bl) for each i G {1, ■ ■ ■ , n}, resulting 
in four possible angles, i.e., wo{i)+wo{iQ)l), WQ{i)+wi{iQ)l), wi{i)+WQ{iQ)l), 
wi{i)+wi{iQ)l). That is, wo{l) and wi{l) can be neighbored with wo{2) and 
Wi(2); Wo{2) and Wi{2) can be neighbored with 1^0(3) and Wi(3); ■ ■ ■ ; Wo{n) 
and Wi{n) can be neighbored with Wo(l) and Wi{l). Hence, there are 0(4n) 
possible angles in total for a given sequence of sub-wedges. We assume the 
angle x+y formed by each pair {x, y) of sub- wedges to be the 'largest' angle in 
a drawing. Then by using the dynamic programming approach of Theorem [2] 
in 0{n) time, we can obtain the smallest angle fn{x,y) in the drawing, and 
hence the aspect ratio for this drawing is {x + y)/ fnix, y). Then optApsRatio 
can be obtained after considering all the 0(4n) possible angles, so the time 
complexity is 0{An x n) = 0{n'^). □ 

Note that the use of dynamic programming allows us to reduce the run- 
ning time of RE2 and RA2 from 0(n^'^) in [o*] to 0{n) and 0(?2^), respectively. 

Theorem 4. The DE2 problem can he solved in 0{n) time. 

Proof. Similar to the proof in Theorem [21 we define 

gi{wtAl),Wt'^{i)) = min x Wt,i2) + Wt'^{2) x wt.^{3) + 

' tje{0,l},2<j<i-l ^ ^ 

■ ■ ■ + Wt'^_^{t - 1) X WtXl)}, 

which can be formulated as the following dynamic programming formula: 

gi{wt,{l),wt'{i)) = min {g,_i{wt,{l) , Wt^ {i - 1)) + Wti_^{i - 1) x wt^{i)} . 

ti-iG{0,l} 

Then, we have 

optSOP = min {gn{wt^{l),Wt'^{n)) + Wt^{l) x Wt'^{n)] . 

Finally, by Equation (|2]), the solution of the DE2 problem can be obtained 
as follows: 



n \ n 

Note that the first and third terms inside the square root of the above equa- 
tion are constants. □ 
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4. Cases C3 and C4 (Unordered Trees with Fixed/Flexible Uneven 
Sub- Wedges) 

In this section, we consider cases C3 and C4 (unordered trees with fixed/flexible 
uneven sub-wedges). For notational convenience, we order all the sub- wedges 
{wo{l),Wi{l), ■ ■ ■ ,WQ{n),wi{n)} in Equation ([T]) in ascending order as 

mi, ma, m„_i, m„, M„, M„_i, Mg, Mi (5) 

where m^ (resp., Mj) is the i-th minimum (resp., maximum) among all. That 
is, Ci = (wtXi) , Wt'.{i)) for i = 1, in Equation ([1]) may be one of the forms 
{mj,mk), {mj,Mk), {Mj,mk), or {Mj,Mk) for some j,k G {1,. ..,«}. For 
convenience, each rrii (resp.. Mi) is said a type-m (resp., type-M) sub- wedge. 

For cases C3 and C4, we consider a bipartite graph G = (V, f/) and a 
function : U f/ — ?■ M in which 

• for case C3, (piV) = {wt-{i) : i = l,---,n}, (j){U) = {wt'Xi) : i = 
1, ■ ■ ■ , n}; for case C4, (P{V) = {Mi, ■ ■ ■ , M„}, 0(f/) = {mi, ■ ■ ■ , m„}; 

• the cost of each edge (f , u) is c(f , u) = 0(i;) + (j){u) for RE, RA and 
DE problems; c(f , m) = 0(f) x 0(m) for SOP problem; the cost of a 
matching N for V x U is c{N) = J2(v u)£N '^('^' '^)- 

Note that, for convenience, each node in ^ U t/ is also denoted by its </> 
function value. 

In case C3 (unordered tree with fixed uneven sub- wedges), for each i = 
1,2, ...,n, sub-wedge WtXi) in V must be adjacent to (matched with) sub- 
wedge Wt'Xi) for some j G {1, 2, n) in U in any solution of our concerned 
problems, and hence the optimal solution must be a perfect matching for 
V xU = {wtXi) : i = 1, ...,n} x {^^/(i) : i = 1, ...,n}. 

In case C4 (unordered tree with flexible uneven sub- wedges) , we have the 
following observation. 

Observation 1. For the RE4, RA4, DE4 or S0P4 problem, there must exist 
an optimal solution in which each type-m sub-wedge is adjacent to (matched 
with) a certain type-M sub-wedge. 

The above observation must hold; otherwise, there must exist k pairs of 
adjacent type-m sub- wedges and k pairs of adjacent type-M sub- wedges for 
some > 1 in the optimal drawing D. But one can easily verify that any 
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(i) Optimal solution 
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55" 43" 40" 30° 25" 
V: Ml Mi Ms m4 mj 
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(ii) Hamiltonian cycle 




/mi —Ml 
(i) Optimal solution 



I h : — N: — 

55" 50" 47" 43" 40° 
V: Ml Ml MpMi Ms 




15" 20" 25" 30" 35° 

(ii) Hamiltonian cycle 



(a) An example for case C3 (b) An example for case C4 

Figure 5: Two examples for expressing optimal solutions as bipartite graphs. 

of our concerned aesthetic criteria of drawing D must be no better than the 
drawing where each of the 2k type-m sub- wedges is altered to be adjacent 
to a certain of the 2k type-M sub- wedges in drawing D (i.e., a drawing in 
Observation [T]). Such an optimal solution in Observation [T] must be a perfect 
matching N ioi V x U = {Mi, M„} x {mi, m„}. 

If Jo denotes the set of the edges corresponding to each pair (wf^(z), Wt'.{i)) 
foriG {!,...,«} (note that (ti'f^(i), Wt'(i)) G^xf/incase C3; (wj. (i), ti;t'(?)) G 
V xVVJV yiUVJU yiVVJU yiU m. case C4), then Jq U forms a Hamiltonian 
cycle for V U . Two examples for the same problem instance but under 
different cases are shown in Figure where the edges in (resp., Jq) are 
represented by dash (resp., solid) lines. As a result, the RE (resp., RA; DE) 
problem is equivalent to finding a matching Nopt ioiVxU such that I^UNopt 
is a Hamiltonian cycle oiVVMJ and the smallest edge cost in Nopt is maximal 
(resp., the ratio of the largest and the smallest edge costs in Nopt is minimal; 
the standard deviation of the edge costs in Nopt is minimal). 

Before showing our results, we introduce some notation as follows. We 
place all the nodes in V (resp., U) on the line y = 1 (resp., ?/ = 0) of 
the xy-plane. Given any matching A^ with two edges ei = (fa, Uh) and 
62 = {vci Ud) inV xU , an exchange on ei and 62 returns a matching A^' such 
that N' = N ® {ei, 62) = {N\ {ci, 62}) U {(fa, ua), {vc, Ub)}. Denote by Cy the 
edge incident to node v in A^. 

Theorem 5. The RES and RE4 problems can be solved in 0{nlogn) time. 

Proof. (Sketch) First consider the RES problem. A careful examination 
reveals that the RES problem and the 2SAL problem are rather similar in 
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Algorithm 2 OptBalloonDrawing-RE3-RE4 
1: construct a bipartite graph VxU — {wti{i) : i = 1, 2, n} x {wt'.{i) : i — 
1, 2, n} for RE3 (resp., V xU ^ {Mi, M2, M„} x {mi, 7712', m„} 
for RE4) 

2: sort the sizes of the sub-wedges in V in nonincreasing order as 

(^2, ■■■,f^n 

3: sort the sizes of the sub- wedges in U in nondecreasing order as 

ai, a2, an 

4: consider a matching N in which ctj is matched with /3j for each i e 
{1,2,. ..,n}. 

5: if /o U A?" is a Hamiltonian cycle for y U C/ then 
6: STOP 
7: end if 

8; order Q = {a^ + /3i+i : i — 1,2, n — 1}, in nonincreasing order 

9: i <(- 
10: repeat 
11: i-(-i + l 

12: if «j and ^j+i belong to different cycles in Iq U N , where aj + ^j+i is 

the i-th maximum in VL then 
13: iV^iV® (e„.,e^.^J 
14: end if 

15: until /o U is a Hamiltonian cycle for y U C/ 
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nature. Hence, Algorithm [2] (a slight modification of the algorithm for the 
2SAL) jsl is sufficient to solve the RES problem in 0(n log n) time. 

The reader is referred to for more details on the proof of the correctness 
of the algorithm. A brief explanation for the correctness is given as follows. 
From [sl, we have the following proposition and property: 

Proposition 1. A matching N determines a solution for RES if IqU N is 
a unique cycle. 

Property 1. Let optAngResl be the optimal solution for RES. Then opt Ang Re si < 
min{/3j + ai,l < i < n}, where V = {/3i, ■ ■ ■ , /3n} ; U = {ai, ■■■,«„}; 
A > ■ ■ ■ > "1 < ■ ■ ■ < ttn- 

See Algorithm [2l If Jq U is a unique cycle at the end of Line 7, then 
Proposition 1 and Property 1 implies optimality; otherwise, Lines 8-15 are 
executed. At each iteration of the loop in Lines 10-15, no matter whether 
N i— N 1^ ,6/3 is executed or not, the cases discussed in [si] can be 
tailored to show that the cost of each matched edge in N is no less than 
optAngResl. Hence, the solution produced by Algorithm |5] must be no less 
than optAngResl. 

The time complexity of the algorithm is explained briefly as follows. It 
is easy to see that Lines 1-8 can be executed in 0{n\ogn) time. At the end 
of Line 7, the nodes of each various cycle are stored in a linked list in 0{n) 
time. Let § be a stack storing the labels top to bottom, in nonincreasing 
order of + A+i. Stack § is used to detect which two cycles we merge next. 
This is done by checking if the endpoints of the edge (oj, A+i), corresponding 
to top element Pi of stack S, belong to different cycles. If they do, the two 
cycles are merged next; otherwise, the element at the top of the stack is 
discarded. Therefore, it takes 0{n) time to detect which cycles to merge. 
The exchanging operation in Line 13 is done in 0(1) time. But also, merging 
two cycles is equivalent to merging two linked lists, which is done in 0(1) 
time as well. As a result, the time complexity of Algorithm [2] is O(nlogn). 

In what follows, we consider the RE4 problem. By Observation [H we 
find an optimal solution for the RE4 problem where each type-m sub-wedge 
is adjacent to a certain type-M sub-wedge, i.e., a perfect matching N for 
V X U = {Ml, M2, Mn} X {mi, 7712, ^n}- By viewing rui (resp.. Mi) as 
ai (resp.. Pi) for each i G {l,...,r7}, the RE4 problem is similar to the RE3 
problem. As a result. Algorithm [2] can also be applied to solving the RE4 
problem in 0(r7logr7) time. □ 
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We now turn our attention to the RA3 and RA4 problems. We consider 
a decision version of the RA3 (resp., RA4) problem: 

The RA3 (resp., RA4) Decision Problem. 

Given a balloon drawing of an unordered tree with fixed (resp., flexible) 
uneven sub-wedges, does there exist a circular permutation a of 
(resp., a circular permutation a of {1, and a sub-wedge assignment t) 

so that the size of each angle is between A and Bl If the answer returns yes, 
then AspRatiOfj t ^ B/A. 

Taking advantage of the analogy between RA3 (RA4) and 2SLW, we are 
able to show: 

Theorem 6. Both the RA3 and RA4 problems are NP-complete. 

Proof. (Sketch) RA3 and 2SLW bear a certain degree of similarity. Recall 
that given a set of n jobs and a range [LB, UB], the 2SLW problem decides 
wether a circular permutation exists such that the workforce requirement 
(i.e., the sum of the workforce requirements for two jobs respectively executed 
at two stations at the same time) for each time period is between LB and U B. 
Given a balloon drawing of an unordered tree with fixed uneven sub-wedges, 
the RA3 decision problem decides whether a circular permutation so that the 
size of each angle (i.e., the sum of two adjacent subwedges respectively from 
two various children) is between A and B. It is obvious that the decision 
version of the RA3 problem can be captured by the 2SLW problem (and vice 
versa) in a straightforward way, hence NP-completeness follows. 

As for the RA4 problem, since the upper bound (i.e., in NP) for the RA4 
problem is easy to show, we show the RA4 problem to be NP-hard by the 
reduction from the 2SLW problem as follows. 

The idea of our proof is to design an RA4 instance so that one cannot 
obtain any better solution by flipping sub-wedges. To this end, from a 2SLW 
instance - a set J = {Ji, J2, Jn} of jobs and two numbers LB, UB where 
Ji = {Wii, Wi2) for each i G {1, we construct a RA4 instance - a set 

of sub-wedges {wo(l), u'i(l), • ■ ■ , WQ{n),'Wi{n)} and two numbers A and B in 
which we let Wmax = max{l^ii, W12, ■ ■ ■ , Wni, Wn2} and p = 2n/ E"=i(W^ji + 
Wj2 + Wmax); woii) = Wa X p and wi(i) = {Wi2 + Wmax) X p for each 
i e {1, n}; A = {LB + Wmax) x p and B ^ {UB + Wmax) x P- 

Now we show that there exists a circular permutation {Js-^, Js2, Js„) 
of J so that the workforce requirement for each time period is between LB 
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and UB a and only if there exist a circular permutation a of {1, and 
a sub- wedge assignment t so that the size of each angle in the RA4 instance 
is between A and B. 

We are given a 2SLW instance with a circular permutation (J^^, J^^, J^^) 
of J so that the workforce requirement for each time period is between 
LB and UB. It turns out that LB < Ws^2 + W^.^ji < UB for each 

i E {l,...,n}. It implies that {LB + Wmax)xp < iW5^2 + Ws,::^,l + Wmax) x P < 

{UB + Wmax) X p for each i G {1, .... n). Consider a = 5 and t = (0, 0, 0) 
in the RA4 instance constructed above. Since wo(c"i) = W^.i x p and 
wi(cri) = {Wai2 + Wmax) X P for cach i e {l,...,n} in the construction, 
thus {LB + Wmax) X P < Wi{ai) + Wo{cri(si) < {UB + Wmax) X P- That is, 
A<9„.<B for each i e {1, ...,n}. 

Conversely, we are given a RA4 instance with a circular permutation a 
of {!,..., n} and a sub- wedge assignment t so that the size of each angle 
in the RA4 instance is between A and B. For any i,j e {l,...,n}, since 

Wl{i) = (Wj2+Wmaa;)xp > Wmax^p > WjiXp = Wo(j)> hcUCC > Wo{j). 

In the RA4 instance, the size of each angle can be Wo(i)+Wo(j), WQ{i)+Wi{j), 
or Wi{i) + for some i, j G {1, ...,n}. For convenience, the angle with 

size wo(i)-l-wo(j) {lesY)-, WQ{i)+wi{j)] wi{i)-\-wi{j)) lox some j G {l,...,n} 
is called a type-00 (resp., 01; 11) angle (note that the order of i and j is not 
crucial here). 

If there exists a type-00 angle in the RA4 instance, then there must exist 
at least one type-11 angle in this instance; otherwise, all the angles are type- 
01 angles. 

In the case when there exists a type-00 angle with size WQ{i) -\- wo{j) so 
that there exists a type-11 angle with size wi{k) + wi{l) for some k, I G 
{1, ■■■,n}, then w.l.o.g., the sub-wedge sequence of the instance is expressed 
as a circular permutation {Si,WQ{i),Wo{j), S2, Wi{k),Wi{l), Ss) where Si - S3 
arc sub-wedge subsequences; the number of sub-wedges in each of Si and 
^3 (resp., 5*2) is odd (resp., even). Let S2 be the reverse of 5*2. Consider 
a new circular permutation {Si,wo{i),wi{k), S2,'WQ{j),wi{l), S3), in which 
the size of each angle is between A and B, because the size of each angle in 
5*3 U Si and S*^ is originally between A and B; A < Wo{i) +wi{k) < B (since 
WQ{i) + wi{k) > wo{i) + wo{j) > A and wo(i) + wi{k) < wi{l) + wi{k) < B); 
similarly, A < wo{j) + wi{l) < B. 

If there still exists a type-00 angle in the new circular permutation, then 
we repeat the above procedure until we obtain a circular permutation S where 
all the angles are type-01 angles. By doing this, the size of each angle in S is 
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between A and B, and the sub- wedge assignment t in the drawing achieved 
by 5 is (0, 0, 0) or (1, 1, 1). In the case oit = (1, 1, 1), we let 5 ^ 5^, 
then t becomes (0, 0, ...0). 

Consider the 2SLW instance (constructed above) corresponding to the 
circular permutation 6. In the 2SLW instance, for each i G {1, ■■■,n}, work- 
force requirement Ws^2 + W^^.eii = {wi{6i) + Wo(5i©i))/p - Wmax- Hence, 
A/p - W^ax < Wsa + < B/p - Wrr,ax, which implies LB < Ws^ + 

W^^.e.i < UB. □ 

We can utilize a technique similar to the reduction from Hamiltonian- 
cirde problem on cubic graphs (HC-CG) to 2SLW ([is']) to establish NP- 
hardness for DE3 and DE4. Hence, we have the following theorem, whose 
proof is given in Appendix because it is too cumbersome and our main result 
for the DE3 and DE4 problems is to design their approximation algorithms. 

Theorem 7. Both the DE3 and DE4 problems are NP-complete. 



5. Approximation Algorithms for Those Intractable Problems 

We have shown RA3 and RA4 to be NP-complete. The results on ap- 
proximation algorithms for those problems are given as follows. 

Theorem 8. Algorithmic is a 2 -approximation algorithm for RA3 and RA4. 

Proof. Let aangResi (rcsp., bangResi and TangResi) bc the minimal angle (resp., 
the maximal angle and the aspect ratio) among the circular permutation 
generated by Algorithm [2j Denote aopt (resp., bopt and Topt) as the maxi- 
mum of the minimal angle (resp., the minimum of the maximal angle and 
the optimal aspect ratio) among any circular permutation. Since bangResi < 
2Mi < 2(x + Ml) < 2b opt where x is the sub- wedge adjacent to Mi in 
the circular permutation with the minimum of the maximal angle, we have 
bangResi < 2bopt- By Theorem |5l we have aangResi = aopt = optAngResl. 

Therefore, r angResl = bangResi/ BangResi < '^bopt/O'opt < 2r opt- D 

Next, we design approximation algorithms for the NP-complete DE prob- 
lems. Here we only consider the approximation algorithms for the S0P4 and 
DE4 problems because the approximation algorithms for the SOPS and DE3 
problems are similar and simpler. Recall that the S0P4 problem is equiva- 
lent to finding a matching Ngpt for bipartite graph V x U , such that c{Nopt) 
is the minimal, where c{N) = u)gjv *^('^) ^ 
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Consider a matching iV^; for bipartite graph VxU in which Mj is matched 
with rrii for each i, i.e., c{Nu) = Yl^=i ^i'^i- Assume that Jq U Nd consists 
of 7] subcycles for 1 < 77 < n, in which we recall that Jq denotes the set 
of the edges corresponding to each pair (wtXi) , Wt'.^i)) for i G {l,---,n}. 
According to matching N^, we have that each subcycle in Iq U N^, contains 
at least one matched edge between Mj and rrii for some i. Let the exchange 
graph x = (^5 -E'x) bipartite graph V x U he a complete graph in which 

• each node in corresponds to a subcycle of Iq U A''^, i.e., = rj; 

• each edge Ci = {u, v) in Ey^ corresponding to two subcycles Cu and 

in Jo U Nd has cost ^/'(ei) = min{ra,bSb,a|(Ma, m^) G (C„,Cj,) U 
(C„, C^) for any a, 6; r^,;, = - M,,, Sb,a = rrib - rria}. (In fact, the 
cost represents the least cost of exchanging edges and mVxU.) 

When ip^Ci) = Vk^iSi^k for some fc, /, we denote /i(ei) = k and z/(ei) = /. Let 
= {V^, Erp^) be a minimum spanning tree over x- With exchange graph x 
and its minimum spanning tree as the input of Algorithm [31 we can show 
that Algorithm [3] is a 2-approximation algorithm for the S0P4 problem. 

Figure [6] gives an example to illustrate how the algorithm works. Fig- 
ure[6]^a) is Iq U Nd where the sohd lines (resp., dash lines) are the edges in 
Jo (resp., in A''^). Figure [6]^b) is its exchange graph x, and we assume that 
Figure M^c) is the minimum spanning tree for x where each edge Cj in 
has weight r^(^ei),uiei)Su{ei),fi{ei)- We illustrate each Si after each modification 
in Line 11 of Algorithm [3] as follows: 
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Algorithm 3 ApproxBalloonDrawing-SOP4 

1: construct the exchange graph x = (^)-^x) V xU 

2: find the minimum spanning tree T-^ = (V^,-Btx) of exchange graph x where 
1^x1=^ 

3: let Si = {M^(e.),m^(e.),M^(e^),m^(e^)} for each edge e,, G Et^ (noticing that if 
i]j{ei) = Tk^is^k for some fc, Z, then /u(ej) = k and v{ei) = I), where each is 
said to correspond to Si (i.e., there are >S'i, S2, • • • , >S'jj-i) 

4: let 5 = {5i,---, Vi} 

5: for each set Sa in S do 

6: for each element x in Sa do 

7: find a set S";, that includes element x but is not considered before 

8: append the elements in set Si, to the end of set Sa (i.e., the dupUcate 

elements are not deleted) 
9: let both edges Cj and Cj correspond to Sa, where edges Cj and ej in 

correspond to Sa and S^, respectively 
10: S^ S\Sb 
11: end for 
12: end for 

13: for each set in S, remove the duplicate elements in each set 

14: order the elements in each set Si, and then denote the new set as S'^ = 

{m'i,m'2, - ■ ■ , m'i,M[, M[_^, ■ ■ ■ , M[} where m[ (resp., M|) is the z-th minimum 

(resp., maximum) in Si; the cardinality of S[ is 21 
15: for each S[ do 

16: M'- is matched with m^+i for j = 1, • • • , Z — 1 
17: M[ is matched with m'^ 
18: end for 

19: output such a matching Napx ioi V xU 
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• Initial: Si = {M2, ma, Mg, mg}, S'2 = {Mi, mi, M7, my}, S'3 = {M5, mg, Mg, 

Si = {M2, m2, Mg, mg}, S's = {M4, m4, Mg, mg}. 

• The elements in 5*4 is appended to the end of Si. 

Si = {M2,m2,Me,mG,M2,m2,Mg,mg}, S2 = {Mi, mi, M7, my}, 5*3 = 
{M5,m5, M8,m8}, S^ = {M4, m4, Mg, mg}. 

• The elements in 5*5 is appended to the end of ^i: 

Si = {M2,m2,M6,m6,M2,m2,Mg,mg,M4,m4,Mg,mg},5'2 = {Mi,mi,M7 
S-i = {M5,m5,M8,m8}. 

Based on the above. Algorithm [3] returns Napx, and Iq U N^px is shown 
in Figure ini^d). In fact. Algorithm [3] provides a 2-aproximation algorithm 
for S0P4. A slight modification also yields a 2-approximation algorithm for 
SOPS. 

Before showing our result, we need the following notation and lemma. 
A permutation vr is a 1-to-l mapping of {1, ...,n} onto itself, which can be 
expressed as: vr = (vr(l), 7r(2), vr(n)) or in compact form in terms oi factors. 
(Note that it is different from the circular permutation used previously.) If 
T^iik) = jk+i for = 1,2, ...,h-l, and Tc{jh) = ji, then (ji,j2, ■■■Jh) is called 
a factor of the permutation vr. A factor with h > 2 is called a nontrivial 
factor. Note that a matching for the bipartite graph V x U constructed 
above can be viewed as a permutation tc : V ^ U. 

Lemma 1. Forn > 2, let X = {xi,X2, ■ ■ ■ ,Xn} (resp., Y = {yi,y2, ■ ■ ■ ,yn}) 
where Xi (resp., yi) is the i-th maximum (resp., minimum) among all. Let 
Q : X ^ Y he a 1-to-l mapping, i.e., a permutation of {1, ■ ■ ■ , n}. // q{X) 
is a permutation consisting of only a nontrivial factor with size n, then 

n n 71—1 

c{.q{,X)) = ^ Xiyg(i) > ^ XiVi + ^ ri^i+iSi+i^i (6) 

i=l i=l i=l 

where ra,b = Xa-Xb, Sc,d = yc-yd for any a, b, c, d. Moreover, z/r^- i+iSj+ij/ - 
Ti^i^iSi^i^i > 1 for each i, j, j' G {1, ■■■ ,n — 1} andj,j'<i, then 

n n—1 

c{g{X)) > ^ Xiyi + ^ ri^i+iSi+i^i + n-2 (7) 

1=1 i=l 

Note that the difference between Equation @ and Inequality (Q is that In- 
equality Q can be applied only when the factor size n is known. 
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Proof. We proceed by induction on the size of q{X). If n = 2, c{g{X)) — 
SLi "^iy^ = ^i?/2 + "^21)1 - XiVi - X2I/2 = ■'^1,252,1 holds. Suppose that the 
required two inequahties hold when n = k. When n = k + 1, 

i(^{i,-,k}\{e-Hk+i)} 

k 

= ^ XiVg'ii) + Xg-i(^k+i)yk+i + Xk+iyg(k+i] - Xg-i(^k+i)ye(k+i) 
1=1 

where g' is a size-/c permutation consisting of a nontrivial factor with size k. 
Then, 

k 

c{g{X)) = + Xk+iVk+i + (a;^-i(fc+i) - Xk+i){yk+i - ye(k+i)) 

1=1 

k 

= ^ Xiyg'ii) + Xk+Wk+i + rg-i(^k+i),k+iSk+i,e(k+i) (8) 

i=l 

For proving Equation (|6]), we replace the first term in Equation (|8]) by 
the inductive hypothesis of Equation (Q, and then obtain: 

k+l k-1 

c{g{X)) > ^ Xiyi + ^ ri,i+iSi+i,i + r^~^k+i),k+iSk+i,gik+i) 
1=1 1=1 

k+l k 
i=l i=l 

since x^-i^k+i^ > Xk and < yt- 

For proving Equation ([7]), we replace the first term in Equation (|8]) by 
the inductive hypothesis of Equation ([7]), and then obtain: 

fc+i fc-i 



c{g{X)) > y^^Xjyi + y^ji^i+iSi+i^i + k-2 + rg-^k+i),k+iSk+i,e(k+i) 

i=l 1=1 
k+l k 

> ^ Xiyi + ri^i+iSi+i,i + k-l 



i=l i=l 



since (0:^-1(^+1) - Xk+i){yk+i - ygik+i)) > - Xk+i)iyk+i - yk) + I by the 
premise of Equation ([7]) (Note that the permutation consists of a nontrivial 
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factor of size n, and hence the case q + = g{k + 1) = k does not occur 
except for n = 2). □ 



Now, we are ready to show our result: 

Theorem 9. There exist 2 -approximation algorithms for SOPS and S0P4, 
which run in 0{n^) time. 

Proof. Recall that given an unordered tree with fixed (resp., flexible) sub- 
wedges, the SOPS (resp., S0P4) problem is to find a circular permutation a 
of {1, ■ ■ ■ , n} (resp., a circular permutation a of {1, ■ ■ ■ , ra} and a sub-wedge 
assignment t) so that the sum of products of adjacent subwedge sizes {SOP^^t) 
is as small as possible. We only consider S0P4; the proof of SOPS is similar 
and simpler. In what follows, we show that Algorithm [3] correctly produces 
the 2-approximation solution for S0P4 in 0(n log n) time. 

From jij, we have c{Nopt) > c^Nd), which is explained briefly as follows. 
From jij, we have that Njy can be transformed from Nopt by a sequence of 
exchanges xi,X2, ■ ■ ■ ,Xn which can be constructed as follows. Let Nk denote 
the matching transformed by the sequence of exchanges xi,X2, ■ ■ ■ ,Xk for 
k < n. We say a node f in ^ is satisfied in A^^ if its adjacent node in N^. is the 
same as its adjacent node in N^. For i = 1,2, ■ ■ ■ ,n, if the sub-wedge Mj is 
satisfied, then null exchange. Otherwise, if the node adjacent to Mj in 

Ni is adjacent to the sub-wedge Mj in Nopt for i j (i.e., Mj is not adjacent 
to rrii in Ni), then let Xi be the exchange between the edges respectively 
incident to Mj and Mj in Ni. Here, by observing each non-null exchange Xi, 
<P{Nopt) - (p{Ni) = TijSj^i > 0. Hence, (t){N,pt) > 0(A^n) = 0(A^d). 

Let 

n 

(9) 

i=l eGEr^ 

We claim that c{Nopt) > clb- Since Jq U N^pt is a Hamiltonian cycle trans- 
formed from IqUNd consisting of rj subcycles, there exist at least rj — 1 times 
of merging subcycles during the transformation (the sequence of exchanges). 
We can view Nopt as a permutation with several factors . There must exist 
a set A of — 1 edges in forming a spanning tree for exchange graph 
X such that each edge in A must correspond to an edge in Nopt which can- 
not be in a trivial factor of permutation Nopt, i-e., it cannot be Mirrii for 
some i. Therefore, by Inequality ([6]) of Lemma [H c{Nopt) > Yll=i ^i^i + 
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Et^ is the edge set of minimum spanning tree of x- 

In what follows, we show the approximation ratio to be 2. Note that 
Napx denotes the matching generated by Algorithm [31 Let S = {j]~lSi and 
i{E>) = U-^e,eET^{fJ^{ei),p{ei)} in Algorithm!! 

n 

'2c{Nopt) > 2ClB > MiTTli + ^ r^(^e),uie)Suie),^lie) 

> J2 {Mf,ie)mf,(^e) + M^(e)m^(e)) + ^ MiTTli + ^ r^(e) 

e&ET^ ie{l,2,---,n}V(S) eGEy^ 

The last inequality above holds since MiTrii for any i G i(S) never presents 
in the first summation term more than twice; otherwise we can find another 
spanning tree with cost strictly less than that of T^. For example, we con- 
sider Figure int^c). Suppose that the cost of edge 63 in is r2,5S5,2, rather 
than r58S8,5, i.e., M2m2 is used three times by ei, 63, and 64 (with costs 
'^2,6'56,25 ''^2,5'S5,25 and ^2,959^25 respectively). We can obtain a contradiction 
by considering a spanning tree T replacing edge 64 by edge C4C5 with cost 
rs 9S9 5, which is less than r2 9S9 2 in general. (The cost of T is less than that 
ofT^O 

Recall that = Ma — Mf, and Sc,d = rric — TUd. Hence, combining the 
first and third terms of the above inequality, we obtain: 

2c(iVopt) > ^ (M^(e)m^(e) + M^(e)m^(e)) + ^ MiTTli 

e&Er^ je{l,2,---,n}V(§) 

i=l i=l ig{l,2,---,n}V(S) 

The above inequality holds due to /i(e) 7^ u{e) for any e G -E^. Since Mgm'j^ > 
M'^gi^Tn'^ in every S'-, we obtain: 

n-i /\si\-i \ 
2c{Nopt) > 5^ 5^ (M;m;.+i) + M|'5,|m; + ^ M^m^ = c(Ar^px) 

*=1 \ j=l / ie{l,2,-,n}V(S) 

In what follows, we explain how the algorithm runs in O(n^) time. 
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In Line 1, the exchange graph can be constructed in O(ri^) time as follows. 
It takes O(n^) time to construct a complete graph x with rj < n nodes in 
which the nodes corresponds rj subcycles in Iq U A^, and the cost of each edge 
is assumed to be infinity. Then, it takes 0{(2)) = 0{n'^) time to compute all 
possible ra,bSb,a = {Ma — Mb) {nib — ^a) for any a,b G {!,■■■ ,n}. Consider 
each TafiSb^a- If Ma and Mi, belong to two different subcycles in Jq U A^, say 
Cu and C^, respectively, and Va^Sb^a < 4^{^i) for their corresponding edge 
Cj = in graph x, then V'(ei) ^ ra,bSb,a- Obviously, after considering all 

possible ra,bSb,a in O(n^) time, graph x is the required exchange graph. 

In Line 2, it is well-known that the minimum spanning tree for graph x 
can be found in 0(n log n) time. Line 3 runs in 0{n) time since each element 
is denoted only once. Line 4 is done in 0{n) time. 

We explain how Lines 5-13 can be done in 0{n) time as follows. Note 
that in Line 3, in addition that each set includes four elements, we record 
that each element knows which set includes it. Hence, in Line 7, any set 
Sb including element x can be found in 0(1) time. Line 8 is done in 0(1) 
time, since each set is a linked list. Note that in Line 7 all the sets that 
includes element x will be considered at the end of Line 12, because in Line 8 
a duplicate element of x is appended to Sa and will be considered again in 
later iteration. Lines 9 and 10 are done in 0(1) time. Therefore, Lines 7-10 
are done in 0(1) time. We observe from Lines 5, 6, 8, 10 that each element 
in S*!, . . . , S'^_i is considered once at the end of Line 12. Since the number of 
elements in 5*1, ■ ■ ■ , Sr,-i is 4(7^ — 1), there are 4(?7 — 1) iterations, each of which 
is done in 0(1) time. Hence, Lines 5-12 are done in 0(4(?7 — 1)) = 0{n) time. 
In Line 13, by scanning each set in S, all duplicate elements are deleted in 
0{n) time. 

Line 14 can be done in 0{n) time, because the ordering of {mi, m2, ■ ■ ■ , m^, 
Mni Mn-i, ■ ■ ■ , Ml} is known. Lines 15-18 are done in 0(n) time, because 
each element is matched only once. □ 

Note that Algorithm [3] is a 2-approximation algorithm for the S0P4 prob- 
lem rather than the DE4 problem because the approximation ratio is incor- 
rect when the minus of the first and third items inside the square root of 
Equation ([2]) is negative. Therefore, we rewrite Equation ([2]) as: 



StdDeVa. 



'j:ti{Mi + m^) , 2J:tl^t'M)'^Ke^{'^i®^) (Y.l=i{M^ + m,) 



TO,; 



2 
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Algorithm 4 ApproxBalloonDrawing-DE4 



The algorithm is almost the same as Algorithm[3]except Lines 15-18 in Algorithm[3] 
is replaced as follows: 

13': for each S[ with \S[\ > 2 (otherwise trivially) do 



14 
15 
16 
17 
18 
19 
20 
21 
22 
23 



let r',=M'-M! and s' , = m' - m' 



an element is said to be available if it is not matched yet 
for each j = 1, • • • , jS^'l — 2 do 

the available maximum is matched with the available second minimum 
else 

the available minimum is matched with the available second maximum 
end if 
end for 

is matched with rn'^.^,^; M^^,^ is matched with m^, where and 
are the remaining elements excluded in the above condition for 
some a, 6 G {1, • • • , \S'j^\ — 1} 
24': end for 



Note that the combination of first and fourth items inside the square root 
of the above equation is the variance of {Mi + mi, M2 + m2, ■ ■ ■ , M„ + ^n}, 
and hence must be positive. Therefore, the DE4 problem is equivalent to 
minimizing the sum of the second and third items, i.e., to minimize 

n n n 

1=1 i=l i=l 

Algorithm m provides an 0(-\/n)- approximation algorithm for DE4. A slight 
modification also yields an 0(-\/?2)-approximation algorithm for DE3. Fig- 
ure [71(a) is an example for Algorithm HI 

Theorem 10. There exist 0{^/n)- approximation algorithms for DE3 and 
DE4, which run in 0{n'^) time. 

Proof. Recall that given an unordered tree with fixed (resp., flexible) sub- 
wedges, the DE3 (DE4) problem is to find a circular permutation a of 
{1, ■ ■ ■ , n} (resp., a circular permutation a of {1, ■ ■ ■ , n} and a sub-wedge 
assignment t) so that the standard deviation of angles {StdDeVa,t) is as small 
as possible. We only concern DE4; the proof of DE3 is similar and simpler. 
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71 = 3 72 = 5 73 = 9 74 = 12 75 = 14 

r'\2 > r 2,3 > r 3,4 < r\s < r 5,6 > r 6,7 > r 7,8 > r 8,9 > r%^w < r'10,11 < r'11,12 < r'12,13 > r'13,14 > r'14,15 

M'l M'2 M'3 M'4 M\ M'b M'l M's M'g M'lo M'n M'12 M'13 M'14 M'15 

\ \ / \ \ \ ^ - " - - / / \ X 

"^^^ >^ ^ ^^-"^ ^--z v"^ 

', ^ , \ , - , , \ , \ , ^ , • , , , . , N , 

mi mi mi m^ ms me ni 7 wi 8 «i 9 m 10 m 11 m 12 m 13 m 14 m 15 



(a) 



r 1,2 r2,3 r3,4 r4,5 r5,6 r6,7 r7,8 r8,9 r9,io r 10,11 r 11,12 r 12,13 r 13,14 r 14,15 




5 2,1 S3,2 5 4,3 5 5,4 S 6,5 S l,b S ij 5 9,8 5 10,9 5 11,10 5 12,11 5 13,12 5 14,13 5 15,14 



(b) 



Figure 7: An example showing how Algorithm |4] works, (a) Certain Ns'. with \S[\ = 15 in 
Napx- (b) Illustration of c(7Vs;) - Y}^^[ M-m'^ induced by (a). 



In what follows, we show that Algorithm H] correctly produces 0{y/n)- 
approximation solution in O(nlogn) time. Let Nopt be the matching for 
V X U witnessing the optimal solution of the DE4 problem, and Napx be 
the matching generated by Algorithm HJ From Theorem [9l c{Nopt) > clb, 
and hence 



n 



n 



c{Nopt) - ^ MiUli I > nicLB-^ MiUli j = n ^ r^(^e),u{e)S,y{e),f,{e) 

v-i /\si\-i \ 

^ (10) 
i=i \ j=i ) 

since /i(e) ^ v{e) for every edge e ^ E^. Observing the matching A^5' for 

each S[ generated by Algorithm H] (e.g., see also Figure [Tl^a)), without lose of 
generality, we assume that 1 < ji < j2 < ■ ■ • < jfc < • • ■ < jfe = I'S'^'I — 1 and 
h is odd such that in 

''M > ?'2,3 > ■ ■ •> ''7i,7i+1 <''7i+lji+2 < • • • < ■Tj^j^ + l 

< r' > r' > • • • > r' 
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M'\ M'lM'iM'AM'^M'b M'lM'iM'iM'AM'sM'i M\ M'lM'i M'a M's M'b M\ M'lM'i M\ M'5 M'e 



I I 



-►X 



^ y \ 



I I 



I I 



rnxmimitrnminK, mumminnmsm^ miimminnmim^ mimirmnnmsm^ 
(al) matching M (a2) matching Ni (a3) matching A^3 (a4) matching Na 



ri.2 



S2A 

C(N2) - c(m) 
= r'i,2i 2.1 

(bl) 



r\.2 r2.i 



S2,l S 3.2 
C(A^3) - C(M) 

= r'l2S'2,l+ r'2.3J 2,1+ r 2,3J 3,2 



r 1,2 r 2,3 r 3.4 



(b2) 




J 2,1 ^ 3.2 i4.3 
c(iV4) - C(M) 



(b3) 



M'l M 2 M 3 M4 Ms Me 
\ V V" / ' 

^ >" \^ I 
\^ \ \ I 

, \ \ / \ I 

m'l m'2 m 3 m 4 m'sm'e 

(a5) matching M 



M'l M 2 M'3 M'4 M'5 M'6 
\ \ I I 



1^1 



m \ m2 mi niA m5 me 

(a6) matching A^e 



ri.2 '■2.3 '■3,4 r4,5 



ri.2 r2.3 '■3.4 ^4,5 rs.6 





S2.1 S3.2 s 4,3 s 5.4 
c{Ns) - c(Ni) 

(b4) 



.S2,l S3_2 S 4.3 S S.4 i 6.5 
C(N6) - c(Ni) 

(b5) 



Figure 8: Illustration of the first several intermediate steps of how to obtain Figure [7I^b) 
from Figure[7lja). For i = 1, • • • , 5, matching iV^+i is obtained by exchanging two edges in 
Ni, as shown from (a«) to (a(i + 1)). (hi) computes c(7Vi+i) — c{Ni), and illustrates the 
relation of the terms used in the cost difference as a bipartite graph, in which each edge 
represents their multiplication relation. 
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Inequality f llip is explained as follows. Since Line 17' in Algorithm H] 
considers the relationship between r'jj^-^ and r^^_]^j_,_2 for j = 1, ■ ■ ■ , IS".] — 2, 
thus, without loss of generality, we use h + 1 numbers (i.e., 1 < ji < j2 < 
■ ■ ■ < jfc < ■ ■ ■ ^ jh = I'S'j'l — 1) to classify all data. Then the data 

is alternately expressed as Inequality f|TT]) . in which r'j^j^-^-i, ...are 

local minimal; r[ 2, r'j^j^+i' ^i4,i4+i' • • • local maximal. 

Then, 

ii-i j2-i 

+•••+( E 

Therefore, 

jl jl i2 j 

= (E E + ( E E '^M+i'^j+ij) ~ ^i2,i2+i'^j2+ij2 

j=i '=j j=ii+i '=ji 

ife+i ifc+i ifc+2 j 

^ ^ (E E '^M + l^ + lj") + ( E E '^M + l^' + lj) ~ ''j/= + 2jfc + 2 + l4fc + 2 + ljfc + 2 

j=ifc j=ifc+i+i '=it+i 

+ ■•• + ( E E<;+i4+u) 

j=3h-l l=j 

Consider Figure El^b) for an example. The above multiplication relationship 
of those r.,. and s.,. for Figure [Tl^a) is given in Figure [Tl^b). Figure |8] shows 
how to transform from Figure [Tl^a) to Figure [Tl^b). 

By Inequality f lTTj) . since r'^ i^^ < t^+^j for j < I < ji or ji < / < j or ■ ■ • 
or j < / < jk+i or jk+i <l< j oi ■■■ oi j <l<jh, we obtain: 

\si\ 

c{Ns{)-J2My^ 
i=i 

j=i i=ii+i 
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Jk+1 3k+2 

+ ■■■ + -3 + lyj^j+is'j+i^j) + ( ^ (j - Jfc+2 + l)r^- j.+i^'+i j) 



-r' s' 

jfc+2>jfc+2 + l jfc+2 + l Jfc+2 




Considering Figure EKb) for an example, c{Ns') - Y.j=i M-m'- < Sr^ 2S2,i + 

2'"2,3'^3,2 + l'"3,4'54,3 + 2r4 5S5 4 + (3+5 — l)r5 gSg 5+4rg 7S7 g+Sr^ gSg^^ 

-'-'"g.io'^'io.g + 2r^Q_;^]^s'^]^ + 3r^j^_;^2'5'i2,ii + (4 + 3 — I)'"'i2,i3'^'i3,i2 + '^^i3M^i4:,i3 + 

l'"l4,15'^15,14- 

By Inequalities ( !T0|) and f lT2|) . we have 

(n \ n 

1=1 / j=l 

In what follows, we explain how the algorithm runs in 0{n'^) time. It 
suffices to explain Lines 13'-24'. Lines 14' and 15' are just notations for the 
proof of correctness, not being executed. In Line 17', j+i and r^+^j^2 can 
be calculated in 0(1) time. Hence, Lines 13'-24' in Algorithm H] runs in 0{n) 
time, because the concerned availability (available maximum, minimum, sec- 
ond maximum, second minimum) is recorded and updated at each iteration 
in 0(1) time (noticing that U and V have been sorted, so has SI); each el- 
ement is recorded as the concerned availability at most 0(1) and matched 
only once. □ 



(12) 



6. Conclusion 

This paper has investigated the tractability of the problems for optimizing 
the angular resolution, the aspect ratio, as well as the standard deviation of 
angles for balloon drawings of ordered or unordered rooted trees with even 
sub-wedges or uneven sub-wedges. It turns out that some of those problems 
are NP-complete while the others can be solved in polynomial time. We also 
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give some approximation algorithms for those intractable problems. A line 
of future work is to investigate the problems of optimizing other aesthetic 
criteria of balloon drawings. 
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Appendix 



• On Proof of Theorem [7] 

Recall that the DE problem is concerned with minimizing the standard devi- 
ation, which involves keeping all the angles as close to each other as possible. 
Such an observation allows us to take advantage of what is known for the 
2SLW problem (which also involves finding a circular permutation to bound 
a measure within given lower and upper bounds) to solve our problems. It 
turns out that, like 2SLW, DE3 and and DE4 are NP-complete. Even though 
DE3, DE4 and 2SLW bear a certain degree of similarity, a direct reduction 
from 2SLW to DE3 or DE4 does not seem obvious. Instead, we are able to 
tailor the technique used for proving NP-hardness of 2SLW to showing DE3 
and DE4 to be NP-hard. To this end, we first briefiy explain the intuitive 



idea behind the NP-hardness proof of 2SLW shown in [13| to set the stage 
for our lower bound proofs. 

The technique utilized in [lij for the NP-hardness proof of 2SLW relies 
on reducing from the Hamiltonian- circle problem on cubic graphs (IIC-CG)0 
(a known NP-complete problem). The reduction is as follows. For a given 
cubic graph G with n nodes, we construct a complete bipartite graph B(y, U) 
consisting of n blocks in the following way. (For convenience, V (resp., U) 
is called the upper (resp., lower) side.) For each node Vi adjacent to fj, Vk, 
vi in cubic graph G, a block Bi of 14 nodes (7 on each side) is associated to 
fj, where the upper side (resp., lower side) contains three f-nodes (resp., u 
nodes) corresponding to fj, f^, and each side has a pair of A-nodes, as 
well as a pair of 6- nodes (as shown in Figure [H]). For the three blocks Bj, Bk, 
and Bi associated with nodes vj, Vk, and vi, respectively, each has a w-node 
corresponding to Vi (because Vi is adjacent to Vj, Vk, and vi). These three 
V- nodes are labelled as Vn, Vi2, and Vi^. In the construction, nodes in V 
and U correspond to those tasks to be performed in stations STl and ST2, 
respectively, in 2SLW. 

As shown in Figure [9l the nodes on the upper and lower sides in Bi from 
the left to the right are associated with the following values 

■ ■ ■ ) -^1,7) — {f^ii i^i ~ ^1 i^i ~ 2, Ki — 2, Ki — 3, Ki — A, Ki — 5), and (13) 
{Bi^i, B,j) = {iK, iK + l,iK + 2, iK + 2,iK + 3, iK + 4,iK + 5),(14) 

respectively, where Kj = (n + 1 — i)K and K is any integer > 7; LB = 



A cubic graph is a graph in which every node has degree three. 
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Figure 9: Illustration of reduction. 





\ 1 




A li- 




ft,? 



Gt+C; Gi+Q 



Ci+a 



Vyl A2i-1 621-1 VH fei V/l A2i 




Vjl Xii-ibii-\ Vki bli vn Xn 




Vjl Xii-ibii-i Vti bli vn Xn 


f 4 • * • f A 




f A f « f « 4 




g A □ f ■ f 


1 \ ' \ ' \ / 
1 \ / \ / \ / 




\i \i \i 1 






1 )( X X 




X X X 1 




' \ / \ A ' 


1 M ; \ / \ 
1 / \ / \ 1 \ 




; \ ; \ 1 \ 1 
/ \ 1 \ 1 \ 1 




1 A \ / \ 1 


i ^ (b d d ^ 




d ^ 6 d 6 d ^ 




1 / \ \ 1 \ 1 

d ^ d 6 d 


vyi Xii-\ bii-i vti bli vn Xn 




vfl Xii-ibii-i vki bli vn Xn 




Vjl Xn-ibii-ivki b'n vn Xn 


(i) 


(ii) 


(iii) 



Figure 10: Three possible transition matchings for B,. 



(n + 1) X - 1 and = (n + 1 ) X + 1 . Each edge \nB{V, U) has weight equal 
to the sum of the values of its end points. 

The instance of 2SLW consists of 7n jobs, in which 2n jobs associated 
with pairs of 6-nodes are Jqi = {(621-1, &2i-i)) (^2i,&2i) ^ 1 < ^ 3n 
jobs associated with t'-nodes are /02 = {(fji) 1*22), ('*^i2, '"ia), {ya^Un) : i — 
1, ...,n}, and 2n jobs associated with pairs of A-nodes are /03 — 
(Aj®i, \i(Q2) : 1 < < Note that Iq = /qi U /02 U /03 is a perfect matching 
for B{V, U), and such a matching is called a cziy matching. 

The crux of the remaining construction is based on the idea of relating 
a permutation of the 7n jobs (J[i], Jp], <^[7n]) in the constructed 2SLW in- 
stance to a perfect matching in B{V, U) in such a way that (Vr[j]2, mod 7n)+i]i), 
1 < i < 7n, are matches. Note that W^[i]2, are the two tasks performed 

by stations STl and S'T2, respectively, simultaneously at a certain time. One 
can easily observe that, because of bounds LB and U B, any matching as a 
solution for 2SLW cannot involve a edge connecting two different blocks, and 
the only edges which can be included in N in each block are the dash hues 
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in Figure M Such a perfect matching N is called a transition matching. If 
IqUN forms a Hamiltonian cycle for B{V, U), then it is called complementary 
Hamiltonian cycle (CHC). 

We use notation (-, ■) (resp., [■, ■]) to indicate an edge of a city matching 
(resp., transition matching). Consider a special transition matching Nd = 
{[Aij,Bij] : i = l,...,n,j = 1,...,7}. IqUNd consists of a master A-subcycle 
Cx = [Ai, A;](A;, A2)...[Ai, A^](A., Aiei)...[A„, A;](A;, Ai), n t;-subcycles Ci for 
i = l,...,n (e.g., Ci = [Vn,uu]{un,vi2)[vi2,ui2]{ui2,vi3)[vis,ui3.]{ui3,vu)), 
and 2n 6-subcycles Cb = 6^) for i = l,...,n. Hence, a CHC for 

B{V, U) is formed by combining the 3n + 1 subcycles. From 13|], in order to 
yield a CHC for B{V, U), there are exactly three possible transaction match- 
ings for Bi as shown in Figure [TUl The design is such that edge (f vi) (resp., 
(fj, Vj) and (fj, ffc)) is in a HC of G if Figure [TOT i) (resp., (ii) and (iii)) is the 
chosen permutation for the constructed 2SLW instance. Following a some- 



what complicated argument, 13|] proved that there exists a Hamiltonian cycle 



(HC) for the cubic graph G if and only if there exists a CHC for B{V, U), and 
such a CHC for B(y, U) in turn suggest a sufficient and necessary condition 
for a solution for 2SLW. 



Proofs of Theorem [71 (Sketch) Now we are ready to show the theorem. We 
only consider the DE4 problem; the DE3 problem is similar and in fact 
simpler. Recall that the DE4 problem is equivalent to finding a balloon 
drawing optimizing optSOP. Consider the following decision problem: 

The DE4 Decision Problem : Given a star graph with fiexible uneven 
angles specified by Equation ([1]) and an integer UB, determine whether 
a drawing (i.e., specified by the permutation a G S and the assignments 
(0 or 1) for (1 < i < n)) exists so that SOP^^t < UB. 

It is obvious that the problem is in NP; it remains to show NP-hardness, 
which is established by a reduction from HC-CG. In spite of the similarity be- 
tween our reduction and the reduction from HC-CG to 2SLW (fill) explained 
earlier, the correctness proof of our reduction is a lot more complicated than 
the latter, as we shall explain in detail shortly. 

In the new setting. Equations f lTS]) and ffT^ become: 

{Ai^i, • • • , Aij) = {K{i), K{i) - 2, K{i) - 3, - 4, K{i) - 6, K{i) - 8, - 9); 
{Bi^i, ■■■ , Bij) = (9m, 9m + 1, 9m + 2, 9m + 3, 9m + 5, 9m + 7, 9m + 9) 
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for i = 1, 2, n where = 9n{2n + 2 — i), n >2, and UB = X]i=i ^I'^i + 
In where Mj (resp., mj) is the i-th maximum (resp., minimum) among the 
14n values. (Note that such a setting satisfies the premise of Inequahty ([7]) 
in Lemma [H and hence can utihze the inequahty.) Hence, we have that: 

Aij > Bij, for any 

Aij > A/c^i and Bi j < B^^i if {i < k) or {i = k and j < I)- 

Note that the above imphes that the j-th upper (resp., lower) node in Bi 
is M^i+j (resp., rriji+j) for i G {l,...,n} and j G {1,...,7}. Define r^^h = 
Mn+a - Mn+b and Sa,b = rrin+a - m^i+b in Bi. Hence, 

n,2S2,i = 2, r2,3S3,2 = 1, ^3,454,3 = 1,^4,555,4 = 4, r^^QSe^^ = 4, rQjSjfi = 2, 

which are often utilized throughout the remaining proof. 

If r2 is a set of transition edges, the sum of the transition edge weights is 
denoted by c{n). UCh = lo^N is a CHC for B{V, U) where Jq (resp., N) is 
the city matching (resp., transition matching) of the CHC and = for i = 
l,...,n (i.e., fiipping sub-wedges is not allowed), then c{N) = J2eeN^i^) ~ 
SOPa,t where c(e) is the weight of the transition edge e. 

Now based on the above setting, we show that there exists a HC for the 
cubic graph G if and only if there exists a CHC Ch = /q U for the instance 
B{V, U) of the DE4 problem such that c{N) < UB. 

Suppose that G has a Hamiltonian cycle Ch- Let Ch = V[i],V[2], ...,f[n]- 
The construction of a solution for S(f/, V") is the same as [l3|, as explained 
in the following. Initiating with there exists a pair (M[2],f[2]) of nodes in 
V xU corresponding to v\2] ^ G because is connected with t>[2]. From (isf . 
we have that C\ is merged with C/, Cj, and respectively in Figure ITU] (i). (ii), 
and (iii). Hence, considering the order of B\^^,B\2], ...,-B[n], in iteration i, by 
choosing the appropriate transition matching, say A/^i,], of Bi from the three 
possible matchings in Figure [TOl A/^^j] merges C[j0i] with the master subcycle 
C\. Besides, since the two 6-subcycles in each Bi also are merged with C\ in 
any matching of Figure [101 we can obtain a complementary cycle traversing 
aU nodes in B{U,V). 

We need to check c{N) < UB. In fact, we show that c{N) = UB as 
follows. It suffices to show that c{Ni) = Ylj=i ^n+jf^n+j + 7 for any i G 
{l,...,n} where Ni is the transition matching for Bi. Denote Ac{Ni) = 
c{Ni) — Ylj=i Mn+jiTLji+j . We can prove that Ac{Ni) = 7 for every matching 
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in Figure dUl Case (i) is shown as follows, and the others are similar: 

Ac{Ni) = M7i+im7i+i + M7i+2?7l7i+3 + M7i+3m7i+2 + M7i+4m7j+5 

7 

+M7j+5m7i+4 + Mji+Qmji+Y + MYi+Tmji+Q - ^ Mji+jmyi+llb) 

= ■^2,353,2 + ^^4,555,4 + re,7S7,6 = 1 + 4 + 2 = 7 

From the above computation, one should notice that if Mj is matched with 
a sub-wedge larger than rrij and Mj+i is matched with a sub-wedge less than 
nij^i for j G {7i + l,7i + 2, 7i + 6}, then Ac(A^j) includes r^j+iSj+ij. 

The converse, i.e., showing the existence of a CHC Ch = IqU N for the 
instance B{V, U) of DE4 with c{N) < UB implies the presence of a HC in 
G, is rather complicated. The key relies on the following three claims. 

(S-1) (Bipartite) There are no transition edges in between any pairs of 

upper (resp., lower) nodes in Ch- 
(S-2) (Block) There are no transition edges in between two blocks in Ch- 
(S-3) (Matching) There is only one of Cj, Ck, and Q merged with the master 
subcycle Cx in each Bi. (Recall that each node Vi is adjacent to Vj, Vk, 
vi in G, and hence the statement implies the presence of a HC in G.) 

For proving the above statements, we need the following claims: 

Claim 1 (see j^) Given two transition matchings N and N' between V and 
U, there exists a sequence of exchanges which transforms N to N'. 

Claim 2 If is a transition matching between V and U and involves two 
edges ei and 62 crossing each other, then c(N) > c(N') for N' = N ^ 
(61,62). 

(Claim 2 can be proved by easily checking c(A^) — c(A^') > 0.) It is very 
important to notice that Claim 2 can be adapted even when Jq U may 
NOT be a CHC. The transition matching where Mj is matched with rrij for 
every j (every transition edge is visually vertical) is denoted by Nd, i.e., 
c(Nd) = YljZi ^j^j- Note that if each edge in A^ is between V and U, we 
can obtain c(A^) > c(Nd) by repeatedly using Claim 2 in the order from the 
leftmost node to the rightmost node of V, similar to the technique in the 
proof of Claim 1 [sj . □ 

Proof of Statement (S-1). Supposing that there exits k > 1 transition 
edges between pairs of upper nodes in Ch, then there must exist k transition 
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edges between pairs of lower nodes in Ch, by Pigeonhole Principle. Select 
one of the upper (resp., lower) transition edges, say ei = (Ma,M5), (resp., 
say 62 = {mp,mq)). Consider N' = {61,62). Then c(iV) - c(iV') = (M^ - 
mp){Mh-mq) > {Mn-runf = I8^{n-lf > 7n. Hence, c{N) > c{N') + 7n. 
By the same technique, we can find N" where each edge in N" is between U 
and V such that c{N) > c{N") + 7kn > c{N") + 7n> c{Nd) + 7n = UB, 
which is impossible. □ 

Proof of Statement (S-2). By Statement (S-1), each edge in the transition 
matching of Ch is between V and U. Suppose there exists at least one tran- 
sition edge between two blocks. Assume there are / blocks, {-B^^, Bk2, i3fej}, 
with transition edges across two blocks. Let kmin — min(A;i, k2, ki). Con- 
sider ci = (Majfrid) is the transition edge between i5fc^.^ and B^. for i e 
{!,...,/}, and kmin 7^ k^. Then there must exist a transition edge connect- 
ing to one of the lower nodes of Bk^.^, say rric, by Pigeonhole Principle, 
and we say the edge 62 = {Mb, rric) where rric and M^, are respectively from 
Bk^^^ and Bk^ for j e {!,...,/} and kj ^ kmin- Note that 62 must cross 
ei because Ma and rUc are in B^^.^, i.e., > M^ and rric < "^d- Be- 
sides, we have Ma > M5 + 9n — 9 and rud > rUc + 9n — 9 because two end 
points of edge belong to different blocks. Consider N' = N ^ {61, 62)- Then 
c{N) - c{N') = {Ma - Mb){md - rric) > (9n - 9)^ > 7n for n > 2. That is, 
c{N) > c{N') + 7n> c{Nd) + 7n, which is a contradiction. □ 

Proof of Statement (S-3). Recall that Iq U Nd in every Bj involves sub- 
cycles Cj, C\, Cfoji-n ^fe' ^b2^■ Ci, C\ from the leftmost to the rightmost. If 
there exists a CHC Ch = IqU N for the instance B{V, U), each 6-subcycle in 
Bi has to be merged with some subcycle in the same Bi by Statements (S-1) 
and (S-2). Ac(iVi) is at least 5 due to the merging of 6-subcycles from the 
following four cases (here it suffice to discuss the merging of 6-subcycles with 
their adjacent subcycles because Ac(A^j) in others cases are larger): 

1. C{,2,_i merged with Cx and Cbj. merged with Ck- ^c{Ni) > r2,3S3,2 + 
'^4,555,4 = 1-^4 = 5 

2. Cb2i_i merged with C\ and C^j. merged with Cf. Ac{Ni) > r2,3S3,2 + 

^•5,6-56,5 = 1+4 = 5 

3. Cb2i_i merged with Ck and Cb^i merged with C^: Ac{Ni) > r^^^s^^^, + 
r4,5S5,4 =1+4 = 5 

4. Cb2i_i merged with Ck and Cb2i merged with Cf. Ac(A^j) > r3^4S4^3 + 

^5,656,5 = 1 + 4 = 5 
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Recall that there are 3n + 1 subcycles in B. Hence we require at least 3n 
times of merging subcycles to ensure these subcycles to be merged as a CHC. 
Since we have discussed that two 6-subcycles have to be merged in each Bi 
(i.e., the total times of merging 6-subcycles are 2n), we require at least n 
more times of merging subcycles to obtain a CHC. In fact, the n times of 
merging subcycles is because each Bi contributes once of merging subcycles. 
As a result. Statement (S-3) is proved if we can show that after merging two 
6-subcycles in each Bi, the third merging subcycles in Bi is to merge one of 
Cj, Ck, and Ci with C\. 

In what follows, we discuss Ac{Ni) when there are exactly h times of 
merging subcycles in Ni'. 

• lih = 2, then Ac{Ni) > 5. 

• If /i = 3 and the transition matching of Bi is one of the matchings in 
Figure [TOl then Ac{Ni) = 7. 

• li h = 3 and the transition matching of Bi is NOT any of the matchings 
in Figure [101 then Ac{Ni) > 7. 

• If 7i = 4, then Ac{N^) > 9. 

• If 7i = 5, then Ac{N,) > 11. 

• If 75: = 6, then Ac{N,) > 13. 

If the above statements on h hold, then Statement (S-3) hold. The reason 
is as follows. Remind that we need 3n times of merging subcycles to be a 
CHC. Therefore, if there exists a transition matching of Bi with h = 2 for 
some i (i.e., there are exactly two times of merging subcycles in Bi), then 
there must exists a Bj for some j with h> A. Then Ac{Ni) + Ac{Nj) > 14, 
which is impossible because this results in the total Ac larger than 7n. □ 

Proof of Statements on h . Note that the transition matching of every 
Bi can be viewed as a permutation of {Mrj+i, M7j+2, ...,Mji^Y} (a mapping 
from V to U), and hence different ordering or different times of merging 
subcycles lead to a permutation with different factors, e.g, the permutation 
for Figure [Hi) is {Mn+i) {Mn+2Mn+3) {Mn+4M7,+5) {Mj^+eMn+r) ■ If we let 
/ = {Mjj^, Mj2, Mji^) be a nontrivial factor of the permutation for A^,, 

then c{Ni) > c{f) > Y.f=n ^krrik + T.i=j{ rk,k+iSk+i,k by Equation in 
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Lemma dJ Here we concern the value X^fctr/j i^k,k+iSk+i,k induced by / (which 
is denoted by 0(/)) because it can be viewed as a lower bound of Ac(iVj). 

If a factor / includes Mj but excludes M^+i, then we say that / has a 
lack at Mj+i. We observe that if the permutation pi for Bi has a lack, then 
we can find a permutation for Bi consisting of the factors without any 
lacks such that Q{p'i) < Q{pi) in which the number of factors of p'^ is the 
same as that of Pi and the size of each factor is also the same. The reason 
is as follows. Assume that pi has a factor / = {..., Mj, Mi, ...) with a lack 
at Mj+i (i.e., I ^ j + 1) and the minimum number appearing in the factor 
is Mq. Let p[ be almost the same as pi except the factor / in pi is modified 
as a factor without any lacks involving Mj+i but excluding Mq in p[. Then 
by Equation in Lemma [H Q{pi) - 6(p-) > (rjj+2'Sj+2,j - ''"jj+iSj+ij - 
rj+ij+2Sj+2,j+i) + Tq-i^qSq^q-i > 2 + 1 > 0. lu the similar way, we can find a 
permutation with factors without any lacks. 

In light of the above, it suffices to consider the permutation for Bi con- 
sisting of the factors without any lacks when discussing the lower bound of 
c{Ni). Thus, in the following, when we say that the permutation for Bi has 
a factor /, this implies that / has no lacks, so Lemma [T] can be applied to /. 

Now we are ready to prove the statements on h. The statement of h = 2 
holds because c{Ni) is increased by at least 5 when two 6-subcycles have to 
be merged in each Bi. As for the statement of h = 6, note that merging 
six subcycles implies a permutation with a factor of size seven. Thus, by 
Equation ([6]) in Lemma [H Ac{Ni) > E5=i = 2 + 1 + 1 + 4 + 

4 + 2 = 14 > 13, as required. Let = J2j=i '^iJ+i-^j+ij = 14 for the 
convenience of the following discussion. As for the statement of /i = 5, 
the permutation involves two nontrivial factors after five times of merging 
subcycles. Note that one of the two factors has size at least four, and hence 
the factor (ji, J4) contributes Y.k=n ^^mu + Yj'k=h ^k,k+iSk+i,k + (4 - 2) 
by Equation ([7j) in Lemma [H Therefore, by Equation in Lemma [H 
Ac(iVi) > ^-r^,:r+iSx+i,x + (4-2) = 16-r^,3;+is^+i,a; for some x G {1, ...,6}. 
(Note that —rx^x+iSx+i,x suggests that and M^+i are in different factors.) 
Since rx,x+iSx+i,x < 4, hence Ac(A^j) > 12 > 11, as required. 

As for the statement of /i = 4, by Equation ([6]), Ac(A^j) > i/j—rx,x+iSx+i,x— 
for some x,y E {1,...,6} and x ^ y. Discuss all possible cases 
of pair {x, y) as follows. Consider one of x,y \s 2 or 3. We assume that 
X = 2, and the other case is similar. Hence, rx,x+iSx+i,x = 1- Since 
Ty^y^iSy^i^y < 4 aud thcrc exists a factor with size at least three in this 
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case, Ac{Ni) >V^-l-4 + (3-2) = 10>9by Equation (jT]), as required. 
The remaining cases are (1, 4), (1, 5), (1, 6), (4, 5), (4, 6), and (5, 6). Consider 
one of x,y is 1 or 6. We assume that x = 1, and the other case is simi- 
lar. Hence rx^x+iSx+i,x = 2. Since Vy^y+iSy+i^y < 4 and there exists a factor 
with size at least four or two factors with size at least three in this case, 
Ac{Ni) > ^Z-- 2 -4 + (4-2) V 2(3-2) = 10 > 9 by Equation ([7]), as required. 
Last, consider {x,y) = (4,5), namely, M4 and M5 (resp., M5 and Mg) are in 
different factors. Hence, M5 cannot be matched with 7714 nor mg, i.e., subcy- 
cle Cb2i cannot be merged with adjacent subcycles Ck, Q. Since merging Cb^i 
with Cb2,_i induces the smallest cost r^^^s^s = 9 in this case, and the other 
two times of merging subcycles must induce cost more than 2, hence Ac(A^j) 
is at least 9. 

As for the two statements of h = 3, by Equation f|T5|) . Ac{Ni) in the case 
when Ni is one of the matchings in Figure [10] is exactly seven, as required. 
Then we consider the case when Ni is not in Figure [10] in the following. 
By Equation ([6]), Ac(iVj) > rx,x+iSx+i,x + ry^y+iSy+i^y + r^^^+is^+i,^ > 5 + 
rz,z+iSz+i,z for some x,y,z & {1, 6} and x ^ y ^ z ^ x since it is necessary 
to merge 6-subcycles, which contributes at least 5. It suffices to consider the 
cases when r^^z+xSz+x.z < 2, which may violate our required. That is, z may 
be 1, 2, 3 or 6. By considering four possible cases of merging 6-subcycles, 
one may easily check that whatever z is, Ac(A^j) must be either larger than 
7 or in Figure [10] □ 
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